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' I present a new approach to the many-body ground state of quantum-HaU systems. The method 

describes the behavior of a two-dimensional electron system at all Landau-level filling factors v, 
^ ' continuously as a function of magnetic field, and enables one to analytically calculate any physical 

Q , value, characterizing the system, at all v. New proposed trial many-body wave functions have a 

■ clear physical meaning and a very large variational freedom, which opens up wide possibilities to 

search for the ground state of the system at all magnetic fields. 

o> ■ 

CN . PACS numbers; 73.40.Hm 



Keywords: Quantum Hall effect (integer and fractional), Strongly correlated electrons, Wigner crystallization. 



I. INTRODUCTION 



The nature of the ground state of a two-dimensional (2D) electron system (ES) in a strong perpendicular magnetic 
field B = (0,0, B) was a subject of intensive theoretical investigations during the past years A great interest 

i to this problem was motivated by the need to explain the integer [|2^ and the fractional |2^] quantum Hall effects, as 
^ ' well as to understand such fundamental physical phenomena as electron-electron correlations at a degenerate lowest 
Landau level and Wigner crystallization of 2D electrons in strong magnetic fields. Recently discovered novel 
intriguing peculiarities of 2D electron transport at higher Landau levels |2^ |3^ provided further impetus to these 
efforts. 



X3 
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O ' The problem of interacting 2D electrons in a magnetic field was discussed in the literature from different points of 
view. First of all, there exists a solution of the classical many- electron problem ||35[-|37|]. If electrons are considered as 
classical point particles, they form, due to the Coulomb repulsion, a regular Wigner-crystal lattice of a triangular, 

l> ■ Ra- a(h+l2/2,V3l2/2) , n,a2 = 2/\/3, (1) 
or a square symmetry, 

00 ■ 

Q ; Rn - a{li,l2), risa^ = 1. (2) 

O . 

' Here li and I2 are integer, and the lattice constant a is related to the average electron density rig. The energy (per 
\ particle) of the triangular electron lattice on a neutralizing positive background with the charge density -l-|e|ns is 
slightly lower than that of the square one, |36y37|] 



^ , ^WCclassrcal ^ _1 ggOSlT, ^WC ,class^cal ^ _1J501325 (3) 

rH ■ 

Q (throughout the paper I will use the i3-independent energy units e^y^nJ/K, where e is the electron charge and n is 
Q . the dielectric constant) . The classical result (^) does not depend on magnetic field. 

Then, there exists the solution of a quantum- mechanical single- electron problem |38-4^]. The single-particle 
Schrodinger equation 

(I use a symmetric gauge) has the solutions Q 

e = = Tt-^c H — ~~~^2^~~^ ' 0<n<cx), — c»<;< -|-oo, (5) 



1 



Here h, c, and m* are the Planck constant, velocity of light, and the electron effective mass, respectively, n and I 
are integer, Uc = \e\B/m*c is the cyclotron frequency, = 21'^ = 2?i/m*Wc, Ib is the magnetic length, and are 
associated Laguerre polynomials. The states with n = and —oo < I <0 belong to the lowest Landau level with the 
energy huJc/2 and the wave functions which I will write in the form 

V'L(r) = 0o,i<o(r) = -^^^ cM-zz*/2), (7) 



where L = —I = 0,1,2, . . . and z — {x + iy)/X is a complex dimensionless coordinate of an electron. 

The quantum-mechanical many- electron problem does not have a complete solution. As the problem is of extreme 
complexity, the main theoretical efforts were concentrated on attempts to guess the ground state many-body wave 
function, at separate points (or in certain intervals) of the magnetic field axis. For a completely filled lowest Landau 
level (the filling factor v — ttX^Us = 1) a trial Hartree-Fock many-body wave function 

*^j = ^det|VL,(r.)|, L,^0,1,...,N-1, (8) 

was proposed in Ref. Q (det^v will denote the determinant of an x matrix), ft has the form of a Slater 
determinant built from the lowest-Landau- level single-particle states (^, and implies that in the A^-electron system 
the states L = Otoi = A^— 1 are occupied by electrons. In the thermodynamic limit N ^ oo the state (|^) is 
characterized by the uniform electron density Q 'T-e(r) — l/vrA^ and the energy per particle |^ (in units e^ ^fnll k) 

CffF = -71-/2 = -1-57080. (9) 

The Hartree-Fock wave function can describe the properties of the system only at the one point {y = 1) of the 
B-axis; otherwise the density of electrons does not coincide with the density of the positive background, and the 
electro-neutrality requirement is violated. 

At a partial Landau-level filling v < \ one could expect an inhomogeneous distribution of the electron density. 
Therefore, first suggestions for the ground state of the 2DES at < 1 were related with the charge-density-wave 
or the Wigner crystal (WC) model The WC many-body wave function |^ 

vI/[^l, = -Ldet|xL=o(r.,R,)|, (10) 



XL(r.,R,) ^ x% = Mr^ - R^.)e— rHBxR,)/0o^ (n) 

is a natural generalization of the classical Wigner-crystal solution: electrons are described by the gaussians ipL=o 
centered at the points Rj of the triangular lattice (0) (here (pQ is the flux quantum). Contrary to the Hartree-Fock 
solution, the wave function (|l^) can describe the properties of the system at all i', continuously as a function of 
magnetic field, as the electro-neutrality requirement is automatically fulfilled due to the relevant choice of the lattice 
constant in Eqs. (|l|)-(|2|). Due to intrinsic triangular symmetry of the solution (p^), one could also expect |^,|| that 
the _B-dependence of the energy of the system contains oscillating functions of 

BSa B 1 

cx — — , (12) 

<Po ns4>Q V 

where S'a is the area of any triangle in the lattice (^. This could lead to certain peculiarities in the energy at i/ = 1/3, 
1/5, and other fractions. Detailed calculations made by Maki and Zotos ||] revealed however no commensuration 
features in the energy (an interesting modification of the model, with less than one electron per elementary cell, gave 
commensuration features at fractional v, but at a slightly higher energy |^,Q). The reason for that is a small overlap 
of the single-particle gaussians '0l=o centered at neighbor lattice points Rj: in order to get the Aharonov-Bohm-type 
oscillations in the energy, electrons should be able to "tunnel" between different lattice points Rj , but this is hampered 
by a negligible overlap of gaussians tpL=Q, centered at different lattice points, at v < 1/3. 
The trial many-body wave functions proposed by Laughlin 0|, 

*n;:s»=n(-^ (13) 
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describe the system at fractional filling factors v = 1/m, where m is odd integer, and are characterized by the uniform 
2D electron density Ug = I/tottA^. Laughlin showed that at v = 1/3 and 1/5 the energy of the liquid state (^3|) is 
lower than the energy of the Wigner crystal. This result was confirmed by exact-diagonalization calculations of the 
energy of a few-electron system |^,^,^ . At = 1 the Laughlin wave function coincides with the Hartree-Fock one (||) 
and has the energy given by ( p|). 

Subsequent theoretical work~p,p^ p7[ in the region < I was devoted to the development and further generalizations 
of the Laughlin theory to other fractional filling factors. It was established that the Laughlin liquid has the energy 
lower than the Wigner crystal at i/ > 1/7; at z/ — the Wigner crystal state has the lower energy. 

In 1996 an interest aroused to the ground state properties at higher Landau levels v ^ I. It was shown that a 
unidirectional charge-density- wave state is likely the ground state of a 2DES at higher Landau levels. This conclusion 
has been confirmed in a number of subsequent publications pO|-p^ and turned out to be in agreement with recent 
experiments 



28-34 



The Laughlin-based approach describes the ground state properties of a 2DES at separate, specific values of the 
Landau-level filling factor. In this paper I present a new method, which describes the physical properties of the 
system at all ly, continuously as a function of magnetic field. I discuss basic principles of the new method, with a short 
demonstration of its advantages, by an example of a few specific calculations. In Section || I discuss possible ways 
of solving the many-body Schrodinger equation and construct an infinite set of trial many-body wave functions for 
the ground state of the system. In Section [II I present an exact analytical method of calculating different physical 
properties of the proposed many-particle states, and show a few examples of calculated characteristics of the system. 
Conclusive remarks are given in Section IV. Mathematical details can be found in Appendices. 



II. TRIAL SOLUTIONS OF THE MANY-BODY SCHRODINGER EQUATION 

A. The problem 

I begin with the many-body Schrodinger equation for a system of N 2D electrons, moving in the plane z = in a 
perpendicular magnetic field B = (0, 0, B), 

7?vl;[^l(ri, r2, . . . , rjv) = i?*[^l(ri, r2, . . . , rjv). (14) 

The Hamiltonian H = K + Vee + Veb + Vbb of the system consists of the kinetic energy term, 

JV 1 ^ / II \ 2 

i=i i=i ^ ' 

the energy of Coulomb interaction of 2D electrons with each other, 

the interaction energy of electrons with the positive background Veb, and of the positive background with itself 
Vbb- Specific expressions for the terms Veb + Vbb depend on the model of the positive background (e.g. randomly 
distributed ionized donors or a uniformly charged jellium disk). The influence of disorder will only briefly be discussed 
in subsequent Sections. In the main part of the paper I will assume that the positive background has a form of a disk 
with a large radius Rb and a uniform 2D charge density -|-|e|7is. Electro- neutrality requires that the radius Rb of the 
disk and the number of electrons in the system N be related by the formula 

nRln, = N. (17) 

Then the energy associated with the positive background is written as 

+ Vbb = -f dr±^^ +Udrf dr' ^'Mr)nbir') ^ 
J "f^k^-i"! 2 7 J K|r-r'| 

where nb{r) — ns9{Rb — r). I will assume that the magnetic field is sufficiently strong (y < 1), so that all electron 
spins are aligned along the z-axis, and will search for an antisymmetric solution of the Schrodinger equation (pT 
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The many-body wave function ^'[^l(ri, r2, . . . , r^r) is a function of A'' vector variables, so that one needs a suitable 
set of basis functions in a 2A^-dimensional space. As in the literature there was an extensive discussion on the 
availability and completeness of different sets of functions in the considered problem p4|-^, I briefly discuss this 
point here. The basis set which I will use in this work consists of the eigenfunctions of the kinetic energy operator 



(ri)0„2/2(i"2) ■■■<!> njv'jv ('^-'v); < Tlj < OO, 



< k < 



(19) 



This set is complete and orthonornial in a 2Af-dimensional space. A properly antisynimetrized subset of these func- 
tions, 



■=det 

/]vT N 



^Li(ri) i/'L2(ri) 



^Liv(ri) 
V'L„(r2) 

'/'L„(rAr) 



< Li < L2 < • . . < < OO, 



(20) 



will be used for expansion of the ground state trial wave functions. The functions (^0|) arc orthonormal in a 2N- 
dimensional space and are the eigenfunctions of the kinetic energy operator K with the eigenenergy NhuJc/'2. 
Now, I discuss two possible approaches to the problem (111). 



B. Perturbative approach 



One could try to solve the problem (^4|) using the perturbative approach. Assume that the Coulomb energy of the 
system (14e + Veb + Vbb) can be considered as a perturbation, and the non-perturbed Hamiltonian is given by the 
kinetic energy operator K. The ground state of the non-perturbed problem is highly degenerate , so that one needs to 
use the degenerate perturbation theory. First, consider the Schrodinger equation ( |l4| ) with N = 1. This is a problem 
of a single 2D electron in the field of a disk-shaped "impurity" with the charge +|e| and the radius Rb = (vrn^)"^/^. 
According to the degenerate perturbation theory the ground-state wave function of the single-electron problem should 
be searched for in the form 



00 

vl/[^=il(r) = VCi^i(r), 



L=0 



(21) 



with arbitrary coefficients Cl- Substitution of the wave function (21) into the Schrodinger equation ( |14| ) leads to a 
secular equation, which allows one to find, in principle, the energy and the wave function of the ground state of the 
one-electron problem. 

Now, consider two electrons moving in the field of a positively charged disk with the radius Rb — (2/7rns)^/^ and 
the total charge -1-2 |e|. An antisymmetric ground-state two-electron wave function is expanded in the orthonormal 
set (|20|) with iV = 2, 



vl/[^=2l(ri,r2) = 



Li=0L2>Li ^ 



V'Li(ri) V'Lsll'l) 

V'Li(r2) iPlA^2) 



(22) 



here Cli,L2 are arbitrary coefficients. Substitution of the wave function (22) into the Schrodinger equation ( |14| ) 
leads to a secular equation, which gives, in principle, the energy and the wave function of the ground state of the 
two-electron problem. 

Similarly, in the problem of N electrons moving in the field of a disk with the radius Rb = (iV/7rns)^/^ Eq. (|l|), 
and the charge density -|-|e|ns, the ground-state antisymmetric wave function should be searched for in the form 



*[^J(ri^ 



r2, 



E 



1 



0<Li<L2<...<Ljv 



/N\ N 



det |V'L, (ri)| 



(23) 



with arbitrary coefficients Cli,L2,....Lm ■ The ground state of the system can be found, in principle, from the secular 
equation resulting from the substitution of the expansion (E3|) into the Schrodinger equation (M). 
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C. Variational approach. Construction of trial many-body wave functions 



Of course, it is absolutely impractical to try to solve the secular equations discussed above at ^ 1. Instead, 
one can use the variational approach. Any trial many-body wave function in the A'^-electron problem can be searched 
for in the form of an arbitrary linear combination ( p^ ) . One of the possible trial wave functions is the Hartree-Fock 
solution (m. It corresponds to the trivial case 



C 



HF 

Li,L2 



N 



(24) 



when the coefficient Cli,L2,...,Ln vanishes at all but one points of the A^-dimensional space of quantum numbers 
{Li, L2, . . . , ijv}- The number of other possible trial wave functions at all v is however infinite, so that a simple 
physical idea would be desirable in order to narrow the number of possible linear combinations in Eq. (|2^) down to 
a reasonable value. 

Such a simple physical idea has been actually used in the literature . Consider the Wigner crystal wave function 
(|lO|). It satisfies the electro-neutrality requirement at all magnetic fields, has a clear physical meaning and corresponds, 
in the limit of very strong magnetic fields — > 0, to the classical picture of distribution of point particles in the 
triangular lattice. One can easily show that the wave function ( |l0[ ) can be explicitly expanded in a series of the type 
(p3|). For the function xo(r,R-) I write 



Xo(r,R) 



-zz' /2-ZZ' /2+z'' Z 



-zz'/2-ZZ'/2 g (^*^)L 
L=0 



L! 



Y^Al{Z)Mv), 



(25) 



L=0 



where Z ^ {X + iY)/\ and 



Al{Z) 

Then, the two-electron function (|l^) is presented as 



Z' 



.-ZZ' 12 



^Zt^ = ^ det 
\/2! 2 



Xo(ri,Ri) xo(ri,R2) 
Xo(r2,Ri) Xo(r2,R2) 



OC OC 



1 



^ ^ Ai,(Zi)Ai,(Z2)^det 

Li=0L2=0 ^ 



AlAZ2) AlAZ2) 



V'Li(r2) '0L2(l"2) 



E E dot 



1 



det 
2! 2 



ii=0 L2>Li 

Quite similarly, for the iV-electron Wigner crystal wave function I obtain 

1 



V'Li(ri) V'L2(i"i) 
'/'Li(r2) ^lA'^2) 



^[^1 - 



E 



0<Li<L2<...<LjN 



detL4L,(Z,' 



det\i}L,{^i) \ ■ 



(26) 



(27) 



(28) 



Eq. ( p^ ) explicitly gives the expansion coefficients Cl^^l2,---,Ln of the Wigner crystal wave function ( pl]| ) in the series 
of orthonormal functions (00) , 



riwc 



detL4L^(Z,) 



(29) 



the coefficients C*)^*^^ are non-zero at almost all points of the iV-dimensional space of quantum numbers 

{Li,L2, . . . , L]si}. 

Here, it is worth discussing the following point. In the Hartree-Fock iV-electron picture, Eq. (|^), electrons occupy 
the states with L = to 



~ N — 1. The last occupied state corresponds to an electron ring with the radius 
''max = A\/ L„iax + 1 = X^/N — Rhy^. At = 1, the radius rmax exactly coincides with the radius of the positive 
background and the electron density ne{r) exponentially decays at (r — Rh) > A. Sometimes, one draws from 
here a conclusion, that at v = I the states with L > Lmax cannot contribute to the expansion ( p3| ) (as the inclusion 
of such states would lead to an essentially non-zero electron density outside the sample) , and hence that the Wigner 



5 



crystal state ( |lO| ) cannot be considered as a trial wave function at v = 1. This conclusion is however erroneous. The 
Wigner-crystal wave function (|l^) is constructed in such a way (|Rj| < Rb for all j, see Figure P, that the electron 
density outside the sample [at (r — Rb) > A] is also exponentially small, but in spite of this, the expansion ( |2^ ) 
contains all the states Li, with < < oo. This remark is valid both for the Wigner crystal state and for all the 
new proposed states discussed below. 

By analogy with the Wigner crystal wave function one can construct an infinite number of trial many-body wave 
functions which satisfy the electro-neutrality requirement at all magnetic fields and have a clear physical meaning. 
Consider the function 



r2, 



■=det 

/M N 



$(ri,Ri) 
$(r2,Ri) 



*(ri,R2) 
*(r2,R2) 



$(ri,Rjv) 
$(r2,Rjv) 



$(rAr,Ri) $(rA,,R2) ... $(rAr,RAr) 



where 



$(r„R,) = (^(r, - R^.)e— r,.(BxR,)/0o^ 
and <^(r) is an arbitrary linear combination of the lowest- Landau- level single particle states, 



oc 



L=0 



(30) 



(31) 



(32) 



It is assumed that all the vectors Rj in Eq. (^0|) are different, and the coefficients satisfy the normalization 
condition J^l I^^P ~ ^- function (30) describes the state, in which electrons are "localized" near or rotate 



around the points Rj uniformly distributed over the 2D plane with the average density Ug. Under this condition 
the electro-neutrality requirement is fulfilled automatically at all magnetic fields. Note that the vectors Rj are free 
parameters of the theory and are not necessarily given by points of the triangular lattice (|l|). They may coincide 
with points of a square lattice (||), of any other type of a lattice, or even be randomly distributed over the 2D plane. 
In general, the configuration of vectors Rj should be chosen in view of minimizing the total energy of the system. 
For example, if effects of disorder are assumed to be important in a given (e.g. low-mobility) sample, the vectors Rj 
may coincide with equilibrium positions of classical point particles randomly distributed over the 2D plane due to 
the influence of remote donors (the influence of ionized donors on equilibrium positions of point particles in a 2DES 
was studied in Ref. ||5^). Considering different configurations of the vectors Rj one can study the interplay between 
effects of disorder and electron-electron interactions in the integer and the fractional quantum Hall effects. 

The wave functions (^0|) imply that electrons localized near different R^-points are in the same single-particle state 
(^). This is also not necessary. One can assume that the single-electron states ip{ri — Rj), localized at different 
points Rj, are different. In this case the functions 4>(ri,Rj) should be supplied by an additional index j. The trial 
many-body wave functions can then be written in a more general form 



vl/[^l(ri,r2,...,r^) = 



1 



det 
Nl N 



$i(ri,Ri) $2(ri,R2) 
$i(r2,Ri) $2(1-2, R2) 



*Ar(ri,RAr) 

$jv(r2,Rw) 



$l(rAr,Ri) $2(rAr,R2) ... $Ar(rAr, Rjy) 



(33) 



where ^j(r^, R^) = fji^^i — R-j) exp[— zTrr^ ■ (B x Rj)/0o] and 

00 



(34) 



L=0 



(the normalization condition J^l l-^i^P = 1 is assumed to be valid for any j). If all the states $j are different, 
the vectors Rj may be the same. Like in (|25|), the functions $j(rj,Rj) can be explicitly expanded in the set of the 
lowest-Landau- level functions iJjl, 



$j(r„R,) 



L=0 

00 

L=0 



-z,-z*/2-ZjZ* /2+z'Zj 



(35) 



L=0 
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where 



oo mm{L,L'} , 



Similar to (|27D~(p8|), one can get from here an expansion of the many-body trial wave function (JS^) in a set of 
orthonormal lowest-Landau- level Slater determinants (EOl), 



oo 

*[^l(ri,r2,...,rw) = V det I Al, ^ det U^^ (r,)| . (37) 



0<Li<L2<...<Lf, 



The trial many-body wave function ( |33D describes a system of particles uniformly distributed over the 2D plane 
with the average density Us- Each of the particles is in the lowest-Landau- level state. If > 1, the states with 
reversed spins, as well as the states from the higher Landau levels may also be included in the expansion (jsj). The 
vectors TLj and the expansion coefficients Bj^ are free parameters of the theory. They are, in general, some functions 
of magnetic field, so that when B varies, the system may undergo transitions from, say, a state with the triangular 
symmetry of the underlying lattice Rj, to a state with the square symmetry of the lattice, or from a state with one 
set of coefficients Bjl, to a state with another set of Bj^. Some examples of such transitions are discussed in Section 



inc. 



It is also worth noting that the more general form of the trial wave function ( p(]| ) or ( p^ ) allows one to overcome 
the problem of negligible overlap of neighbor single-particle states in the Wigner-crystal solution (^0|) at small i^. As 
one saw from the discussion in Section |, the underlying triangular symmetry of the Wigner-crystal state could lead 
to commensurate energies at certain values of but did not because of the small overlap of the neighbor gaussians in 
strong magnetic fields. Using other linear combinations of the functions V'l in Eqs. ( p^ ) or ( p^ ) (with Bjl 7^ ^£,0)1 
one can get a large overlap of the neighbors keeping the symmetry of the underlying lattice Rj. This may lead and 
does lead to oscillating magnetic-field dependencies of the electron-electron interaction energy of the system, with 
deep minima in certain intervals of v. 

From this point on I will consider a few examples of possible trial wave functions, assuming that (i) effects of disorder 
are negligible, so that the vectors Rj form a lattice; (ii) the lattice of the vectors Rj has either the triangular (|l|) or 
the square (|^) symmetry; (iii) the coefficients Bjj^ in the expansion (^) do not depend on j, i.e. all single-particle 
states centered at different lattice points are the same, and (iv) the expansion ( ^4[) consists of the only one term, 
BjL' = Thus, the trial many-body wave functions to be considered in this paper are 

vI/fU_L.det|xL(r„R,)|, (38) 

where the functions xl are defined by Eq. ([TTI). The functions ^'^^ are supplied by an additional subscript L which 
labels different variational solutions of the problem. They describe a system of electrons, localized at {L = 0) or 
rotating around {L > 0) the points Rj of the triangular or the square lattice. The state 5*^0 is the Wigner crystal 

state The functions XLi^^i, Rj) with different Rj are not orthogonal to each other, so that the trial functions 5*^' 
take into account electron-electron correlations in an approximation far beyond the Hartrce-Fock one. 

It would certainly be meaningless to discuss the complicated many-body wave functions in the form ( ^8|) or (^3|), if 
there were no possibility to calculate the energy or other physical properties of these states. Fortunately, there exists 
an exact analytical method of calculating expectation values of practically any physical quantity, including the density, 
pair-correlation function, energy, et cetera, in the states (|38|), (Eot), or (^3|). I present this method in the next Section. 



III. PROPERTIES OF THE MANY-BODY STATES 

In order to calculate expectation values of different physical quantities in the states one needs to know the 
norm (^'^'|^'^'), and the matrix elements of arbitrary single-particle and two-particle operators. This can be done 
with the help of general results formulated in Appendix 

All the specific calculations below have been performed in a finite-size Wigner-cluster configuration, with the radius 
of the positively charged disk Rb, determined by Eq. (pT|), and N lattice points (|l|) or (||), placed within a circle of the 
radius Rh, Figure |l|. Configurations with N = 7, 19, 37, . . . and N — 5, 13, 29, . . . lattice points have been considered 
in the triangular and the square lattice, respectively. 
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A. The norm 



According to (A4), the norm is determined by the formula 

(*[^l|vI/[^l)=det|5"L 



(39) 



where S^^ is the N x N matrix of overlap integrals with the elements 



S. 



{Xa^\Xa,) = exp 



Here are the Laguerre polynomials, the complex variables are related to vectors = {Xi, Y^), 

and -qij = ii- £,j. 

The norm (\E'^'|\E'^') exponentially tends to zero with growing N, 



(40) 



(41) 



(42) 



where a depends on the overlap integrals, see Figure |^. This property does not lead to any problem in calculating 
the expectation value of a physical quantity O, as in the formula O = («'^^l|0|4'^^V(*L^'l*L^'> both the nominator 
and the denominator contain similar integrals, both tend to zero in the limit N oo, while the ratio O tends to a 
constant. 



B. Single- particle operators 

Matrix elements of single-particle operators 



N 

6i=^oi(r,) (43) 
1=1 



can be found with the help of Eqs. (A6) and (|A7]). In order to calculate them one needs to know the matrix 



Sj^ (Section [II A) and single-particle matrix elements {Xai\oi{ra)\Xai) of a specific operator 6i. Due to a simple 



analytical form of the functions xl, the single-particle matrix elements {Xai\'^i{''Ca)\Xaj) can be analytically calculated 
for practically any physical quantity of interest. A few examples are presented below. 
Let 6i be the kinetic energy operator. For the single-particle matrix element I get 

(XaAlx^,) = ^J drxl{r,Il^) (p + x XL(r,R,) = (44) 



so that immediately 



(45) 



the kinetic energy per particle in the state is Tiujc/'^- 

A less trivial example of a single-particle operator is the electron density 



i(r)-^<5(r-r.)=5]nGe^°-, (46) 



G 



G are reciprocal lattice vectors. It can be calculated either directly or via its Fourier transform. In the former case 
one needs the matrix elements 

^''(r) = {XaMr - ra)\x^^ = x*(r, R.)x(r, R,), (47) 
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in the latter case - Fourier components 



S"(q)^(Xa.|e^''-"HXa,) 



exp 



+ (e*+*772)fe+z7/2) 



Ll [{m,-il/2M^+ij*/2)], 



(48) 



where ^ — {qx + iqy)^- Together with (A6) and (A7), Eqs. ( ^ ) and (^ ) give exact analytic expressions for the 
electron density and its Fourier components in the states . 

Figure |^ exhibits coordinate dependencies of the normalized electron density, rie(r)/ns, in the states ^ l=o and 
^L=2, with a triangular lattice of vectors Rj, at a few values of the Landau level filling factor v. The Fourier 
components of the electron density have been calculated for a large number of configurations (Figure left panel) 
with up to ~ 1000 (they are with a good accuracy proportional to l/ViV), and then extrapolated to — > oo. 
Figure ^ thus shows the electron density in the thermodynamic limit. An interesting finding is that, in spite of the 
underlying periodic symmetry of the lattice Rj , the density of electrons Ue (r) = 1/ ttA^ is uniform in the limit v = 1 
of a completely filled lowest Landau level. This is due to a strong overlap of the neighbor single-particle states at 
V — \, and seems to be valid at all L, at least for a triangular lattice of points Rj. I have calculated the square-root 
deviation of the electron density from the average value. 



1/2 



A, 



El 



LL 



(49) 



as a function of magnetic field B for a few lowest L-states in the system with a triangular lattice of points Rj , Figure 
^. The functions Al^o{B) have a complicated oscillating dependence on the magnetic field, due to the overlap of 
electron rings rotating around the lattice points, and all the curves Ai(_B) tend to zero at v 1. 

Similarly, one can calculate expectation values of other single-particle operators in the states L at any magnetic field. 
Results of a more complete analysis of the single-particle properties of the states will be reported in subsequent 
publications. 



C. Two-particle operators 

Matrix elements of two-particle operators 



N N 



^2 = E E °2(r„r,) (50) 



can be calculated with the help of Eqs. ( |Aq ) and ( A9). In order to calculate them one needs to know, apart from the 
matrix 5^^, the two-particle matrix elements (A9) of a specific operator 02. For instance, the matrix elements of the 
pair-correlation function can be easily calculated. I will focus below on the most important two-particle operator in 
the problem - the electron-electron interaction energy Vee , Eq. ( p^ . 
According to Eq. ( |A8| ) one gets 

N N N N 

iKei^^i) = iEEEE(-i)^^^'^'^'«g-(^-^>g^^(j--o^4" detj5^^i™:f4,., (51) 
1=1 j=i fe=i 1=1 

where the (A^ — 2) x {N — 2) matrix Sj^^2\coimnn:jij i formed from the matrix S^^ , Eq. (^), by deleting the rows 
k and the columns j, I. The matrix elements V^^^ are determined by 



I. 



ijkl 

Vl^k = / / dr,dnix^,rxaj ,^ ^^^^^ {xLTxL (52) 



have the following symmetry properties 

^ijkl ^ klij 



^ijkl — ^klij^ "^ijkl — {^jilkY ^ (53) 
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J\f — 2 ' "'^ \ column^i,j 



and are calculated in Appendix With the help of (|53|), Eq. ( |5l| ) can be presented in the form 

TV N 

N N N 

+ E E E (-iy+'sgn(*-j)sgn(z-fc)Re 

W TV TV N 

+ oE E E E (-ir+^+'+'sgn(z-A;)sgn(j-ORc 



ll'jiifc - ^ikji) ^^^\'^N~2\colunin^i,k 



^ijkl |Ow_9L„;... 



Af-2 



A^ — 2 I column^ j .1 



(54) 



where the first, second and the third Unes correspond to the two-site (direct and exchange), three-site and four- 
site contributions to the total energy of the system, respectively. Dividing this result by the norm (|39|), I get the 
expectation value of the energy of electron-electron interaction in the state , 



N N 



N N N N 

E E y^^~lT. E 



1 ^ " ^ " ~rLL 



Hartree terra 



N 



EE E (-lF+'sgn(z-j)sgn(z-fc)Re 2(y,f/,-yJj^,)a:i 



N N N N 



+ JE E E E (-l)^+■'■+'+'sgn(^-fc)sgn(J-ORe[^,^,V;■/] 



(55) 



where the first (underbraced) term is the Hartree energy, see Appendix and the rest - the exchange-correlation 
contribution. The factors in (pq) are determined by 



A + \ qLL \row^i,k 

detAr-2 \-JN-2\colurnn^j,l 



(56) 



Electron-electron interaction energy (Vee)^' grows as N^^'^ at iV — > oo (the energy of a charged disk). In order to 
compensate this growth, I subtract the energy of the positive background 



Vbb = 



e\/W:8N N 



K 3 V TT 



(57) 



and divide the result by N (to get the energy per particle) . Thus calculated estimate of the electron-electron interaction 
energy 



ef(iV) 



{Vee - Vbb)^L^ 

N 



(58) 



differs from the cohesive energy of the system by the term {Veb + '2.Vbb)^L^ / N . With the help of Eq. ( p5| ) the energy 
e|f (A^) can be further rewritten as 



e|'=(iV) = '^'■"■^'^(TV) + e^^"™''''(iV), 



(59) 



where the Hartree term 



Hartree 



(60) 



includes the energy of the positive background, the exchange-correlation energy 
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N N 



^ex — corr 



(N) 



'■=1 J=l(7^») 




AT N 



N 



E E E (-1) 

TV W N 



sgn(i - j>gn(z - fc)Re 2(y^/^ - Vmog 



N 



91V VvE E E E (-i: 
■'=1 fc=i(#» j) i=i(#ij,fc) 



sgn(i - fc)sgn(j - ^Re 



y~'LL ik 



(61) 



consists of two-site, three-site and four-site correlations, $ — B/ns(j)o ~ 1/i^ is the magnetic flux per particle in the 
units of flux quantum, the matrix elements Vjj^i are related with Vjj^i by the formula 



* tjki 



(62) 



(see Appendix and all the energies are measured in units yfrTsj k. 

The Hartree contribution (|6^) can be calculated exactly in the thermodynamic limit N ^ oo 



(Appendix |^), 



Hartree 



lim e 

N~*oo 



Hartree 
L 



(N) 



'^V^ E 7 




(63) 



Here G are reciprocal lattice vectors, the function !Fl is defined by Eq. (C5), and the integral is calculated in Eq. 
(|C7|). The exchange-correlation contribution is calculated, in practice, with a finite number N of lattice points. It 
converges when N exceeds some number Nq, where A^o = ^oi^, v) varies from a few tens to a few hundreds, dependent 
on the overlap of the neighbor single particle wave functions xl- 
Equations (p8|), (63) and (|6l|), together with 



i)-(B7) and (C7), provide exact closed-form analytical expressions 
for the energy of electron-electron interaction in the many-body quantum- mechanical states V^l, at all values of the 
Landau-level filling factor v, continuously as a function of magnetic field. However complicated, these expressions give 
a principal possibility to calculate the energy of the system in the thermodynamic limit with any desirable accuracy. 

Figure || exhibits the magnetic field dependencies of the Hartree contribution to the energy of the states "^l, 
calculated with Eq. ( |63| ) (thick curves), for a few lowest L, for the triangular (Q) and the square (H) lattice of points 
Rj- [thin curves at L = 1 show, for comparison, the results of a finite- A^ calculations with the formula (|60|)]. As 
expected, the Hartree energy increases with L, and in the limit of the very strong magnetic field all the curves tend 
to the classical- Wigner-crystal line. Note that the influence of the quantum-mechanical effects, even in the Hartree 
approximation, is much stronger than the difference between the energies of the classical triangular and the square 
Wigner-crystal lattices, Eq. The quantum-mechanical giomiA state may thus be realized by the states (|30| ) or (|3^ ) 
with different types of the underlying lattice. Which type of the lattice corresponds to the true quantum-mechanical 
ground state may also depend on magnetic field. Therefore, different types of the lattice should be tried in the search 
of the true ground state of the system at different v. 

Figure || exhibits the magnetic field dependencies of the exchange-correlation energy (|6l| ) of the states 5* l , calculated 
with a finite number N — 7, 19, and 37 of triangular lattice points, for the states with L = 0, 1 and 3. One 
sees that the energy g^^^-<:°rr^j^^ negative and grows with N in its absolute value. Figure]^ shows how the energy 



ex—corr 



(N) converges with A^. As seen from this Figure, up to a few hundreds particles should be taken into account 
in order to get a reasonable accuracy of the final result at the considered values of L and i'. If the overlap of the 
neighbor single-particle states is not negligible, the magnitude of the energy e 

^Hartree 



ex — eorr 



(N) can be comparable with 

(A) (note that the number of Hartree terms grows as A^, while that of the three-site and four-site terms - 
as A'^ and A'^'* respectively). The magnetic field dependence of the exchange-correlation energy has a complicated 
oscillating character, with pronounced minima, resulting from the interference effects of electron waves and the overlap 
of electron-ring states XL(r,Rj). Qualitatively this can be seen in Figure ^d, where the insets show the overlap of 
electron rings at a few characteristic points of the magnetic field axis. The points A and C, for instance, correspond 
to the touch of the nearest-neighbor and the next-nearest-neighbor rings, respectively. In the point B any three 
nearest- neighbor rings intersect in one point. The depth of the minima, and their position on the i?-axis, depend on 
the index L and the symmetry of the lattice. 

Figures p|-pl] show the total energy of electron-electron interaction eff (A^) = e|f'"'*'''^^ -I- eff"'^"^^ {N) in the states 

The Hartree contribution (the uppermost curve) has been 



with L = to 3, as a function of $ = B/ns4>o 
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calculated exactly in the thermodynamic hmit, the exchange-correlation contribution - with a finite number N of 
lattice points. The total energy e|f also oscillates with magnetic field, due to the physical effects discussed above. As 
a consequence, the ground state of the system may be realized by the states with different L (and different symmetry 
of the underlying lattice), in different intervals of magnetic field. 

Results of a more complete analysis of the energy and other two-particle properties of the states '^l will be reported 
in subsequent publications. 



IV. CONCLUSIVE REMARKS 



I have presented basic principles of a new method in the many-body theory of quantum-Hall systems. Only a small 
part of specific results (for the density and the electron-electron interaction energy of the states \E'l) have been shown 
with illustrative purposes. A complete study of all properties of the system in different quantum-mechanical states, 
and all possible consequences from this method will require substantially more time, place and efforts. 

The new approach presented in this paper, involves two essential constituents: (i) a physically clear and transparent 
way to construct an infinite number of trial many-body wave functions, for all magnetic fields B (Section and (ii) 
an exact analytical method of calculating any physical property of the proposed many-body states (Section |lll[ ). With 
the help of this method, one can get comprehensive information on all physical properties of the ground, as well as 
excited, states of the system at different values of magnetic field. 

The new approach has a huge variational freedom and opens up new fields of research in the many-body theory of 
quantum-Hall and other many-body systems: 

1. It can be used for a systematic search of the ground state of a 2DES at all i^, both in the fractional and in the 
integer quantum-Hall regimes (the generalization to the case of non-spin-polarized systems and to the region of higher 
Landau levels has been briefiy discussed in Section ilC). 

2. It makes it possible to study the interplay between effects of disorder and electron-electron correlations in both 
the integer and the fractional quantum Hall effects (Section ilC). 

3. It offers an opportunity to separately analyze the Hartree, exchange, and correlation contributions to the total 
energy of the system at all magnetic fields (Section [II C). This may find applications in the density functional theory. 

4. The new approach can be generalized to and applied in other many-body problems, such as, for instance, many- 
electrons quantum dots, quasi-one-dimensional electron systems, 2D electrons in inhomogeneous external (electric and 
magnetic) fields, metal-insulator transitions in 2DES, and other. 
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APPENDIX A: MANY-BODY MATRIX ELEMENTS 



Let ip{r,a) and 0(r,/3) be arbitrary functions dependent on the 2D coordinate r and arbitrary (maybe, vector) 
parameters a and (3 respectively. Let 



1 



*Ar(ri,r2, . . .,rN;ai,a2, ■ ■ ■,aN) = ^= det \^p{ri,aj)\ 



m N 



(Al) 



and 



<I>Ar(ri,r2, . . . ,r7v;/3i,/32, ■ • ■ ,^w) = ^L= det |0(ri, ^j)| 



Nl N 



(A2) 
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are N x N Slater determinants formed from these functions, and ai, i = 1,2, . . . , N , are two sets of N different 
parameters (a^ 7^ a^, f3i ^ f3j). The following theorems can be straightforwardly proved. 
Theorem 1. The scalar product 

(^atI^at) = J *^(ri,r2, . . . ,rAr;ai,a2, ■ • ■ ,aAr)^'Ar(ri,r2, . . . ,rAr;/3i,^2, • ■ • ,^Ar)dridr2 • • -drjv (A3) 
is determined by the formula 

(*w|$A,) = det^AT = det\S{a,,Pj)\ , (A4) 
where Sn is the N x N matrix with the elements 

S{a,,l3,) = = / r{r,a^)(b{r,f3J)dr. (A5) 



Theorem 2. Let Oi be an arbitrary single-particle operator (^3|). The matrix element {'i> n\Oi\^n) is determined 
by the formula 

N N 

(vl/^|6i|cf>^)= ^^(-1)^+^01 (a.,/3,0detJ5^_J^J„^„^^, (A6) 
i=i j=i 



where the {N — 1) x {N — 1) matrix ^jv-il^ohlmnT^j formed from the matrix Sjy, Eq. ( A5), by deleting the row i 
and the column j, and 

oi{a^,Pj) = j ip*{r,a.{)di{Y)(t){Y,l3j)dY. (A7) 

Theorem 3. Let O2 be an arbitrary two-particle operator ([50|). The matrix element (^Ar|02|<I'Ar) is determined 
by the formula 



N N N N 

(vt^lO^I*^) = E E E E(-l)^'''''''''«gn(* - ^)^g^0' - 0o2(a., afc, /3„ A) drt \Sn-2\7oZtn^,,i^ (A8) 
1=1 j=i fc=i 1=1 



where the {N — 2) x {N — 2) matrix '5Ar_2|™;™^'^^j ; is formed from the matrix Sn, Eq- (A5), by deleting the rows 



i, k and the columns j, I, the function sgn(ri) is defined so that sgn(O) = 0, and sgn(ri) = ±1 for n > (n < 0), and 
02{ai,ak,l3j,Pi) = / V*(i"i,ai)V'*(r2,a/c)o2(ri,r2)(?!)(ri,/3j)0(r2,/3i)c^ridr2. (A9) 



APPENDIX B: COULOMB MATRIX ELEMENTS 

In this Section I present results of calculation of the Coulomb matrix elements 



V^ii = / / dr,dr,(x^,)*x^^,^- -lixLTxti, (Bl) 

where, for generality, I assume that the indexes L and L' are not necessarily equal to each other. In all the formulas 
below S^i are defined by Eq. (^), r/,y — £,i ~ and F^^^z) is the n-th derivative of the fimction 

F(z) = e-'Ioiz), (B2) 

where /q is the modified Bessel function. 

There exist five different types of the integrals V^^[ , dependent on whether one or two pairs of subscripts equal to 
each other, or all the indexes ijkl are different: 

1. Direct Coulomb two-site matrix elements V,^/:; , 
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2 kX L\L'\2^+^' ^ -p)\ V 8 
^0 V 8 ) ^\{2p - ^)\ 



Here and in Eq. ( |B4[) below, = ?7??*/4, and 77 = 77^-. 
2. Exchange two-site matrix elements K^,^ , 



2kA L!L'!2^+^' ^p!(L-p)!V8 



(2p)!(-2)' ^ (2g)!(-2V 

Z^ fc!(2p-i-fc)!(2g-j-/fc)!C'= ^ ' 



k=0 

3. Direct three-site matrix elements Vjj'i^i , 



V 2 2 ; 



V^H - \/ 2 ^ ^ '''L!r'!(-2)'^+^' '''''' 



inin{L+p,L +9} 



(L+p)!(i' + g)! 

X 



E 



X 



,^0 H(L+p-r)!(L' + (7-r)! V??^ir,4^ 

iL-p)l fju_\ ' ^^p" (^'-g)! (^ + P-0! 

^ z!(L-p-z)! V2r^,J _ ^ _ ^■)! 1^2,7*J Z. + p - r - /c)!2fe 

L'+q-r . . , inin{i+fcj+/} 

Z^ /llT' -I- « - -r - Z^ 



/!(L' + g-r-0!2' ^ sl{i + k - sy.{j + I - s)l \ 7^uV*k 



l=Q ^ ^ ' s=0 

X 



4. Exchange three-site matrix elements V^^^ , 

^ E ^i(r^.^-^^^ E )j(-2r E ki{L-ky. {'2^, 



X 



Z^;!(L'-0! 1^ 27y*J Z. s!(^ + fc-s)!(J+^-s)! Vr?.,r?.fe 
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X 



Z^^p\[L-k-p)\\ A^^J ^^q\{L' -l-q)\\ A7^*J 



^ [L - I + p)\{L' - ] + q)\ ( 4 



r=0 



r\{L - i + p - r)\{L' - j + q - r)\ \riijri*,. 



p(L'-j+q+L-i+p-r) I VijVik j ^ ^gg^ 



5. Four-site matrix elements V^^i , 



T.ll{L'-iy\ 4h j ^„ p! (fc - p) ! V '7,7^/*,. 



inin{p,Z/ — 



„xRp-,-i} - 'Z)!(i' - ^ - 'Z)!(J + ^ - P + 9)! V ^fc J ^0 H(z - r)! V^/ry; 
xi.C<+w;-P-r)(^^,^*j4). (B7) 



ik 



APPENDIX C: HARTREE APPROXIMATION 

The Hartree approximation works well when the overlap of the neighbor single-particle wave functions is negligible. 
In this approximation the many-body wave function corresponding to ( |3^ ) is written as 

*L^Lrtree(l"l,r2,...,rA,) = X11X22 ■ ■ ■ Xn N ■ (CI) 

The functions x are normahzed. Hence, {'i'^L,Hartree\'^^L, Hartree) = 1' ^^'^ 

N N 

in ,Hartree .Hartree/ 

i=l j = l 



where T^^^ is given by (p^). 

The density of electrons in the state ^'L.ffartree is 



N 

L, Hartree f 

N—^oo ■ 

1=1 



= (r)= hm ^|V'L(r-R.)P^5]n^«''^'-^(G)e^°-, (C3) 
where G are reciprocal lattice vectors, 

N 

L,Hartree^Q^ = Us / lim V " R.) Pe-^°-"-dr 

|V'L(r)pe-^°--dr = n.^F^ ((GA)V4) , (C4) 
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and 

=e-"LL(x). (C5) 
The formula for the Hartree energy can be transformed as follows 

i=l j=l{=ii) J J \ a b\ 

= H I {"e ■'''"■*''''(^<^)"e '""""^^^(r.) - n,{r,)n,{n) - N\Mra)\'\Mn)n (C6) 

where nb(r) = nsO{Rb — r) is the density of the positive background. Substituting here the Fourier expansion from 
( |C3|) I get Eq. (^) from the main text of the paper. The integral there can be analytically calculated, 
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FIG. 1. Examples of finite-size Wigner-lattice configurations which were used in calculations: left plots - triangular lattice 
with N = 19, 37, and 187 lattice points, right plots - square lattice with N = 13, 29, and 149 lattice points. Circles show the 
boundary of the uniform positively charged jellium disks. 
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FIG. 2. Logarithm of the norm log{^'i|>I'_L) as a function of the number of lattice points involved in calculations, for a few 
different L and v, for a triangular and a square symmetry of the lattice. 
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FIG. 4. Square-root deviation Al of the electron density from the average value Us, Eq. (p9[), as a function of the inverse 
Landau-level filling factor "I> = = B/nacfio, for a few lowest ^'^-states in a system with a triangular lattice of points Rj. 
The results are obtained by extrapolation to the thermodynamic limit A'^ — > oo as described in the text. 
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FIG. 5. Hartree contribution to the energy of the states as a function of the inverse Landau-level filling factor 
calculated in the thermodynamic limit for a few lowest L, for a triangular (a) and a square (6) configurations of lattice points 
Rj. The energies of the classical Wigner crystal ^ are shown by thick dashed lines. Thin dotted curves (L = 1) exhibit the 
Hartree energies, calculated with a finite number N > 2000 of lattice points using Eq. (|60|). These curves illustrate that the 
Hartree energy converges very slowly with A'': to get a few-percent accuracy, one needs to take into account a few thousands 
particles. 
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FIG. 6. Exchange-correlation contribution to the energy of the states "Pl as a function of the inverse Landau-Ievel filling 
factor calculated with a finite number A'^ of (triangular) lattice points N = 7, 19, and 37, for (o) the states L = and L = 1, 
and (6) the state L = 3. Insets in Figure b qualitatively show configurations of electron rings at magnetic fields corresponding 
to the energy minima in a few characteristic points A, B and C of the l/i'-axis. 




-2 I ■ ' ■ ' ■ ' ■ 1 

50 100 150 200 

number of lattice points, N 

FIG. 7. Calculated exchange-correlation energy per particle e2"^~°°'"''(iV), in units e^y^nT/zt, at z/ = 1, as a function of the 
number of lattice points N in the state with L = to 3. These curves illustrate the convergency of the exchange-correlation 
contribution to the energy of the states ^l- to get a few-percent accuracy one needs to take into account from a few tens to a 
few hundreds particles, dependent on L and u, i.e. dependent on the overlap of neighbor single-particle states. 
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FIG. 8. Calculated electron-electron interaction energy per particle e^^iN) = ^H^rtree _|_ ^^^-'^"'■'"(tv), in units e^i/n7/K, as a 
function of B/ns4>o = 1/v in the state 5'l=o (the Wigner crystal, Ref. [H]). The Hartree contribution (the uppermost curve) 
is calculated exactly in the thermodynamic limit, the exchange-correlation energy - with a finite number N = 7, 19, 37, ... of 
the (triangular) lattice points. 
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FIG. 9. The same as in Figure 0, but for the state 
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FIG. 10. The same as in Figure ^, but for the state "i! l=2 
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FIG. 11. The same as in Figure H, but for the state '^l^s- 
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